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Abstract 


It is a simple consequence of the Cafarelli-Kohn-Nirenberg theory that 
every possible singularity in a thin Haussdorff-measurable set of a Leray- 
Hopf solution of the incompressible Navier Stokes equation is on the tip of 
a small open cone, where the solution is smooth. Using global regularity 
results for weak Hopf-Leray solutions this potential singularity can be 
analyzed by investigation of the asymptotic behavior at infinite time of a 
solution of a related initial-boundary value problem posed in transformed 
coordinates on a cylinder. It follows that the velocity components and 
their spatial first order derivatives are left continuous at each potential 
singular point. Next to some new consequences such as as global regularity 
of the Leray Hopf solution after finite time (for L 2 -data) many known 
results can be recovered with this method succinctly, for example the 
result that H 1 -regularity implies global existence and smoothness, and 
that local regularity follows from L 3 -regularity of the data. 


1 Some simple observations about Leray-Hopf 
solutions in the light of the CKN-theorem 

Hopf told us essentially that for L 2 -data h = (hi, ■ ■ ■ , h n ) and positive viscosity 
v > 0 the Navier-Stokes equation initial value problem with periodic boundary 
conditions, i.e., the problem 





(1) 


< div v = 0, 


, v (0, •) = h, 


has a global weak solution for the velocity components Vi, 1 < z < n, where 


£ L°° (R+, L 2 (T*)) n L 2 oc (R+, H 1 (TP)) . 


( 2 ) 


The pressure p is determined by the velocity components via the Leray projec¬ 
tion. Here, the symbol T n denotes the torus of dimension n and R + denotes the 
set of positive real numbers. The arguments in [2], or in [3], tell us in addition 
that 

i) there exists an open dense subset / C (0, oo) such that tul (R+ \ I) = 0, 
where mi denotes the Lebesgue measure, and such that 


Vi eC°° (I x T”), 1 < i < n, 


(3) 
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where C°° denotes the space of smooth functions. 


ii) 

S := {(t,x) G R+ x T n |Vi C°° at (t,x) for some 1 < i < n} (4) 
has vanishing one dimensional Hausdorff measure. 

iii) for any T > 0 and for dimension n = 3 a weak Leray-Hopf solution 
Vi, 1 < * < 3 satisfies 

Vi&Li ([0,T],L 4 (T 3 )) (5) 


for all 1 < i < n. 

Proofs of these results can also be found in standard text books such as volume 
III of Taylor’s book on partial differential equations. The open dense set I may 
be represented by an union of open intervals, i.e., 

I = L)j£jUj, (6) 

where J is an index set and Uj C R+ are open intervals in the field of real 
numbers R equipped with the standard topology. Next if t s G R+ \I is a time 
of a time slice {t s } x T” related to a possible singularity, then there exists an 
index j € J and an open interval Uj = (to,t s ) C I. According to [2j we know 

Vi G C°° (Uj x TP) for 1 < i < n, (7) 


while (O implies 

Vi G ([t 0 ,t s \,L A (T 3 )) , and v t G L§ (P„L 4 (T 3 )) (8) 

especially. The gap to a proof of uniqueness is filled if for any given finite horizon 
T > 0 we can prove that a Leray-Hopf solution Vi , 1 < i < n satisfies 

i) Vi G L°° (R +,L 2 (¥")) n L 2 ([0, T], H 1 (T*)); 

ii) Vi GI 8 ([0,T],I 4 (T 3 )). 

From the perspective of weak function spaces it seems difficult to close the gap 
between item ii) and (151) . while the condition in item i) is quite close to what 
is known by Hopf’s result. In this context recall that in the case of dimension 
n = 3 the standard (mollified equation) arguments for Hopf’s result in © 
use compact sequences with limit Vi G L 2 ([0, T], 7J 1_e ) for any e > 0 and 
T > 0. Indeed, these standard arguments tell us that a family of solutions 
iP, l<*<n, fc>lof mollified equations build a compact sequence in 
L 2 ([0, T], i7 1_e °) for arbitrarily small eo > 0 such that (passing to a subsequence 
denoted again by (v £ i k )k> 1, 1 < i < n if necessary) the set 


Ite [0, T]|V 1 < * < n |u I £l (t,.)| ffl _ eo +^2\v e i k (t,.)-v e i k "(t, Oljyi-eo <°o \ 

[ /c>2 J 

(9) 

is dense in [0,T]. This leads to the open dense time set where the weak Leray- 
Hopf solution is smooth, such that any possible singularity after any finite time is 
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at the endpoint of an open interval (to, to + A to ) of smoothness (for some A to > 
0). More precisely, for a zero sequence (e*,)*, the weak limit of (v e i k )k> 1 , 1 < i < 
n satisfies 


Vi = lim v e , k 
fefoo 


£ L 2 ([to, to + A t0 ], H 1 -^) n C°° ((to, to + A t0 ) x T"). 


( 10 ) 


Hence possible singularities at (t a ,x s ) for positive time t s > 0 are at the larger 
endpoint of such an open time interval (as is well-known). 

Next for given t s > 0 assume that (t a ,x s ) £ S' is a singular point of a 
Leray-Hopf solution. For given j £ J consider the corresponding open cone 

K (ts,Xs) {(t,x)\t £ Uj = (to, is) k\x- x a \ < t s - t} , (11) 

where |.| denotes the Euclidean distance. For this j and for all 1 < i < n we 
have Vi £ C°° ^ ta ’ Xa) ') and 


for all (t,x) £ Kf‘ ,x,) 


\vi(t,x)\ < yj- 


c 

t)^\x — a; s | A 


( 12 ) 


for some finite constant c £ and for some parameters /x, A which satisfy (case 
n = 3) 

0<m<|o<A<| (13) 

Concerning first order spatial derivatives we have for all 1 < i < n Vi £ 
L 2 ([ 0, T], iJ 1_e ) for small e > 0 (by the standard mollified equation arguments), 
hence 


for all (t, x) £ Kj t ‘ > ’ Xa ' > 


u i>fc (t,x)| < 


(t s 


c 

t)^° |a: — x s | A ° 


(14) 


for some finite constant c £ R+, and for some parameters Mo? Aq which satisfy 


1 3 

0 < /xo < ^; 0 < Aq < - + e for small e > 0. 


(15) 


Here, K^ s ’ Xs ' > denotes the closure of the open cone K^"" Xe ' > 


i.e., 


K (ts,x B ) g [to, t s ] & |x - x a \ < t g - t} . (16) 

Next we introduce a technique to push the singulatity orders m and A of the 
velocity components Vi,l < i < n from | and | respectively to zero and, 
similarly, the singularity orders m 0 and Ao of the first order spatial derivatives 
Vij, 1 < i,j < n from A and ^ respectively to zero. This implies that the 
velocity components and the first order spatial derivatives of the velocity com¬ 
ponents of a Leray-Hopf solution are left continuous at each potential singular 
point (t,x) £ S. Since there may be singularities in the slice {to} x T”, we 
consider the restricted domain 

■= {(*,*)11 £ (ti,tg) & ti £ Uj \x-x s \ <t s -t}. (17) 

Then for all parameters which satisfy (fTSll the functions : Kj t ^ Xs ' > — > R, 1 < 
i < n defined by 

tq ,M (t, x) = ( t s - ty\x - x s \ x Vi(t, x) (18) 
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and their first order spatial derivatives uf’j 1 ,1 < i,j < n are bounded functions 
on the closed cone K^f ,Xs \ where 

3,t l 1 

u* £ C°° , 1 < i < n. 


Next we define classes of singularities. 

Definition 1.1. A function / : [0,T] x ffi™ — > R is said to have an isolated 
singularity from the left, i.e., with respect to increasing time, if there is an open 
cone Kj t £’ x *' > ( defined above) such that the restriction / < t s ,x s ) of / to Kj t ^’ Xa ' > 

J ! ^ 1 

is smooth and has a singularity of order (A, p) for positive real numbers A, p at 
the point (t s ,x s ). Here, a smooth real-valued function f K (t B ,x s ) £ C°° (k^^ 3 ^ 
is said to have a singularity of order (A, p) for positive real numbers A , p at the 
point (t s ,x s ), if for some t\ £ Uj and cone K there is a finite constant 
c > 0 such that for all (t,x) £ we have 


\f{t,x)\ < 


\t — t s \v\x — a; s | A ’ 


(19) 


and if there are no X < A and p' < /i and a finite constant d such that 


\f(t,x)\ < 


1 1 - t s |^'|a; - x s | A ' ‘ 


( 20 ) 


The singularity analysis below shows that the velocity components u* them¬ 
selves and their first order and second order spatial derivatives have no singu¬ 
larities from the left (of the type just described). Furthermore, we shall show 
that there is a uniform upper bound (from the left) at each time section t s x P. 
This implies that a Leray-Hopf solution has a left-continuous extension. We 
shall combine this with local time regularity arguments in order to obtain global 
regularity results after any finite time for the Hopf-Leray solution. As a first 
consequence of the CKN-theory and singularity analysis we have 


Theorem 1.2. For a given time horizon T > 0 let Vi, 1 < i < n be a weak 
Leray-Hopf solution of the Navier Stokes equation, where £ L 2 ([0, T], H l ) 
for all 1 < i < n. Then after any finite time this solution Vi, 1 < i < n has 
no singularities of order (p, o,Ao) from the left for any 0 < Xo,po- In other 
words, the Leray-Hopf solution is left-continuous an the time interval (0,T] for 
arbitrary given T. Moreover, the assumption can be weakened assuming that Vi £ 
L 2 ([0, T], iJ 1 ” 6 ) for any e > 0, and after finite time there is still no singularity 
of order (/xo,Ao) from the left such that the relations in \15\) are satisfied with 
po > 0 or with Aq > 0. 


The theorem is proved in the next section. Next we consider the main 
idea of singularity analysis considered in this article, and state some further 
consequences. For each possible singular point (t s ,x s ) of a given weak Hopf- 
Leray solution Vi, 1 < i < n there is an open cone 3 where the Leray- 

Hopf solution component functions are smooth. We consider the coordinate 
transformation c : K^f' Xa) -> x (ts ’ Xs) , where 

JXl f i ’ 

r, z) =(^,—L_), p £ (0,1.5). (21) 


4 







Here p is a parameter. For the analysis of this paper it is sufficient to assume that 
p = 1 (such that z[ t i s ’ Xa ' > becomes a cylinder). This choice has the advantage 
that the transformed equation has neither (weakly) degenerate (for p < 1) 
or singular (for p > 1) second order coefficients. This implies that except in 
the case p = 1 the singularity analysis needs to be extended by an analysis 
of fundamental solution of heat equations with degenerate or weakly singular 
coefficients. Therefore, for simplicity, we shall stick to the case p = 1 in this 
paper. Note that Z t v is a cylinder of infinite height. Given a global weak 

Hopf-Leray solution and for a potential singularity at time t s we choose ti £ 
Uj = ( tj,t s ) for some j £ J and consider that global weak Hopf-Leray solution 
locally in transformed coordinates on the cylinder 

Zg°’ xa) = [t in ,o o)xQ, (22) 


where [fj n , oo) := 


t *_) ti , oo) . For 1 < i < n and all (r, z) £ zj: t i a ’ Xa ' > define 


Wi(r,z) = Vi(t,x). 


(23) 


Remark 1.3. Note that we have a family of comparison functions Wi = w^ s ’ Xb ' > 
here (for each possible singularity ( t s ,x s ) we construct one). If it is clear from 
the context that we refer to a given (t s: x s ) we drop this superscript for the sake 
of simplicity of notation. 

We have 


dzj 1 

Vi j = Wi i — — = Wi i -r - 

,J ,J dxj ,J (t s - t)P 




(ts~typ- 


(24) 


The function Wi , 1 < * < n is more regular with respect to time. Nevertheless, 
as Wi(r,z) = Vi(t,x ) we can transfer information obtained for Wi , 1 < i < n 
to Vi, 1 < i < n. From the perspective of CKN theory this is an advantage: 
assuming H l regularity we have to push the regularity from Ts ( Uj,L 4 (T 3 )) 
to L 8 (Uj,L 4 (T 3 )) in order to obtain uniqueness, and if we can do this for 
Wi,l < i < n, then we can do it for the velocity components Vi, 1 < i < n 
of the correspeonding Leray-Hopf solution itself. Moreover, if we can weaken 
the assumption of H l -regularity to IL 1_e -regularity in this context, then we 
have uniqueness after any finite time in the situation of Hopf’s theorem. More 
precisely, if for all possible singularities (t s ,x s ) with t s > 0 and for a cylinder 
zf,t( Xs) = [Un,t s \ xllwe have 

Wi £ L°° ([ t in ,t s \,L 2 (fl)) n L 2 ([ti,i s ], H 1 (H)) n L 8 ([t in ,t s ], L 4 (H)) , 

then after small time to > 0 we have 

Vi £ L°° ([t 0 ,T],L 2 (T n )) nl 2 ([0,T],i? 1 (T n )) n L 8 ([0, T ], L 4 (TP)) . 


Note that for strong data, say for data in H m r\C rn , m > 2 the latter statement 
implies uniqueness and regularity for all time because we have a local time 
contraction result for such strong function spaces at initial time t = 0. Therefore, 
even from the the perspective of weak function spaces the study of the more 
regular function Wi, 1 < i < n has some advantages. For this reason we analyze 
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the behavior of this function on a cylinder which corresponds to a cone in original 
coordinates. Next, we derive initial-boundary value problems for each cone with 
an assumed singularity at the tip of a related cone. First we observe that the 
incompressibility condition is conserved on the time interval [U n ,t s ), since 


n 

0 = div v = Vi t i 

i= 1 





(25) 


For the time derivative we have (recall that p = 1) 


dr x j 

Vi , t = w itT -^ + 2 ^ w <3 _ t y > 

3= 1 


(26) 


where for t £ [0, t s ) we have 

dr _ d t _ 1 


t _ _ tx—t+t _ _ ts 


dt dt (ts-t) (ts—t) ' ( t B -t ) 2 (ts—t) 2 (ts — ty 


> 0. 


(27) 


Hence the function wt , 1 < i < n is determined by initial-boundary value 
problem 

^ - M2«'E"=1 w i,j,3 + Ml E”=l + Zj) Wij = -piV iP w , 

div w = 0, 


K (t s ,x s ) , 


(28) 


d°K yz . 

3,t\ 


where Q s z^“' Xa ' > denotes the spatial boundary of the cylinder Z^‘ ,Xs \ and 


M2 = M2(r) = £(t a - t) 2 2p , Mi = Mi( T ) = - ^) 2 p , (29) 

and where t = t(r) denotes the value at r of the inverse r —> t(r) of r = r(t). 
Note that we used 

Zj = (t s — t) p Xj. (30) 

Furthermore, p w denotes the pressure in transformed coordinates. Note that 
is bounded for our choice p = 1 in this article and becomes weakly singular for 
p e (1,1.5). 

Remark 1.4. We shall later show that we have a regularity transfer from w — 
i, 1 < i < n to Vi, 1 < i < n, i.e., full regularity of the comparison functions 
Wi, 1 < i < n (for each cone one) implies full regularity of the original velocity 
function Vi, 1 < i < n. Note that for the sake of notational simplicity we drop 

^(*s,a:s) > x s) 

the reference to the cone in general, i.e., we have Wi = w i 1 = v i 41 for 

(t X ) 

each cone K tx , and where for all 1 < i < n v i 1 denotes the restriction 
of Vi to the cone K^‘ ,Xe \ For refined regularity investigations we may use a 
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variation of initial-Neumann-boundary value problems of the form 


TFT - M2^ E"=i w ij,j + Mi E"=i ( w j + z i) w i,j = -Mi Vip w , 


div w = 0 





(31) 


where d u denotes the normal spatial derivative. Indeed in case p = 1 we gain 
regularity of one order for the estimation of the boundary terms in the estimation 
of Wi if we use Neumann boundary condition (factor (t s — t) p ). 

Remark 1.5. The term ’weakly singular’ means -roughly- ’integrable’. Further¬ 
more, note that lim 7 -.|. 00 p\{r) = 0, and lim T -|- 00 pq(t) = 0, and these coefficients 
are bounded in any case. We shall analyze the role of the coefficients at time 
t = t s in case p G (1,1.5) elsewhere, where we reconsider the Levy expansion of 
the fundamental solution in this context. 

The problem described in (12811 is a initial-boundary value problem. The 
initial time may be abbreviated by 


tin '■= --—, where t s > 0. 


(32) 


In the case p = 1 and p 2 = j~ is a constant. Hence, the fundamental solution 
G p of the transformed heat equation 


(33) 


q,r - P'lV&q = 0 


is an explicitly known Gaussian type function (no expansion is needed). Apply¬ 
ing the divergence operator to the first equation in (12811 we get 



(34) 


as the incompressibility condition transfers to the function Wi, 1 < i < n, and 
where we recall that we add the superscript w to the pressure in order to indicate 
that the pressure is considered in transformed coordinates. The Leray projection 
form of (BHD is determined via 


n 



Hence the Leray projection form is obtained from (l28l) if we replace the first 
dynamical equation by 



(36) 
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Here, L\ n denotes the Leray projection operator in case of the n-torus (corre¬ 
sponding to the ith derivative of the pressure). This operator can be determined 
explicitly by Fourier transformation, which we shall do below in the singularity 
analysis. 

Hence, in the domain zj: t i s ' Xs ' > (resp. the cutoff domain Z := [U n ,T\ x H) 
we have a the classical representation 

m(r, z ) = f yeQ Wi(t in , 2/)G()(t tm, z y)dy 
- f z (ts,* s ) (e;=i Ml (Wj + Zj)wi^ (s, y)G^(r - s,z - y)d?/ds 

. . x ( 37 ) 

+ /xi(s)L^.„ (^Ey m =i (s, y)G^(r -s,z- y)dyds 

+ /ti f d s z «s,*s) wf(s,y)GZ(T,z-,s,y)dyds. 

1 1 > T 

Here, for the boundary term we use (IdOD below. In the following the abbreviation 
TV* refers to the sum of the Burgers term and the Leray projection term (cf. (TUI) 
below). Then the boundary terms are determined by 

w? Z (t,z) = -2vi\ SK (T,x{z)) +2 j n w i (t in ,y)G l i,(T,z-,t in ,y)dy + ( 2TV* * GO) (r,z) 
EfcLi /tL In ( - 2 ^| Sif (<b *(0) + 2 f n Wi{t in , y)G>}(a, £; t in ,y)dy 
- (2TV* * G(j) (a, £)) x G0’ fc (r, z, a, 

(38) 

where is defined in (1321) above, S'x denotes the spatial boundary of the 
cone, z —> x(z) denotes the transformation from the spatial boundary of the 
cylinder d s Z^ ta,XB ^ to the spatial boundary of the cone Kj f ^ Xs \ and the series 
G(j ,fc , k > 2 is defined recursively by 


GO’ 1 := GO, 


G0’ fc+1 (r,u;s,u) := J Q G^(T,u;a,w)G^(a,w,s,v)dS y da. 
Furthermore, note that 


(G | QS Z ( l s < x s) Vi | (t, X) , 


Sk ' 


(39) 


(40) 


and 


(2TV* * GO) (t, z) = 

~ 2 f z (u,*,) (E;=iMiK +^)w*,j) {s,y)G»{T - s,z - y)dyds ( 41 ) 

+2f z (H,*.) Mi(s)L^. n (E",m=i Kj%,m)) (s, 2 /)G 0 (r -s,z- y)dyds 

Here, recall that GO is the fundamental solution of the equation G0 r — 
y^^AGO = 0. Furthermore, the coefficient p 2 is bounded for p = 1, but for 


p £ (0,1) or p £ (1,1.5) we have a ’degenerate’ or ’weakly singular’ coefficients 
P 2 respectively. 

Using such classical representations of the transformed comparison function 
Wi, 1 < * < n, Theorem P is proved in section 2. For each argument (t s ,x s ) 
of a possible singularity of a velocity component functions we find an upper 
bound and we can show that for each time t s we can find an upper bound for a 
family of comparison functions which do not depend on x s . As a consequence 
the superior limites of the velocity component functions have a lower and upper 
bound although we cannot construct this upper bound by the methods in [2]. On 
this abstract level we can argue: if there is no upper bound of the modulus of a 
velocity component function, then we can find a point in the Hausdorff set with 
a (A, /^-singularity which satisfies 0 < po < 0 < Ao < §.. However Theorem 

P tells us that this is not possible. We have even more. The CKN-theory tells 
us that for the set S of singularities we have 

H\S) = \imH^(S)=0 (42) 

where for 1 < p < n 

H pS = |£ ( diamB i) P : S c u 7= i B i and diam(Sj) < <5 j . (43) 

Let S ta := {(t,x) £ Sjt = t s }. We have S ts C S and 

H n (S ts ) < H\S ts ) < H\S) = 0 (44) 

n / 2 

and it is well-known that for = Aj^-T (n/2 + 1) 

7 n H n is Lebegues measure. (45) 

Hence sharpening Theorem 11.21 for derivatives of the functions Wi, 1 < i < n 
with uniform upper bounds and transfer of the result to the original velocity 
components Vi , 1 < i < n leads to the possibility of regular extension of these 
functions to the time section {t s } x T” for all finite t s > 0. Note that for 
multivariate spatial derivatives of finite order a = (aq,--- ,a n ) we have for 

M>i 

\D*Vi(t,.)\ <c a \D°Wi(T,.) l \ pla | , (46) 

where |.| denotes the spatial supremum norm, i.e., |/| := sup x6R n |/(a;)| for 
a function / : R" —> R, and |a| := Hence, if we estimate multi¬ 

variate spatial derivatives of the velocity components of a Leray-Hopf solution 
using local solution representations as in (ED it becomes increasingly diffcult to 
get a regular upper bound as we aim at upper bounds for higher order deriva¬ 
tives. Nevertheless we can transfer regularity results for the comparison function 
Wi , 1 < i < n to the corresponding original Hopf-Leray solution using spatial 
Lipschitz continuity of the Euler-Leray data function applied to regular data 
(which is the application of the Leray projection term operator to given data). 
We may interpret the initial value problem for 17 ,1 < i < n as an initialbound¬ 
ary value problem, where we consider the data on the boundary of a cylinder 
to be given by the solution. We then have a similar representation as in (1371) 
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above. More precisely, for some constants 0 < Ci < C 2 < 1 and in an interval 
\t\,t s ) the original velocity function Vi, 1 < i < n has a local time representation 
on a cylinder C1 C2 := {(t,x)\t £ [t\,t s ) & 2 ; £ [ 01 , 02 ]} of the form 


Vi ( t , x) = /| y ! (ti , y)ez?i ci c2 } Vi (h , 2 /)G„ (i, x-y)dy 

- fz« (Z?=i v j v ij) ( s > y)G v {t -s,x- y)dyds 

+ / z - (E? m =i ( v m,jV jtm )) (s,y)G„(t- s,x~y)dyds 

"ti.cj_.C2 \ / 


(47) 


+ /f t , /osz” v? Z (s,y)G v (t,x;s,y)dyds, 

1 t1.e1.c2 

where G„ is in case /z = 1, and where the boundary term has an analogous 
definition as in the case of the comparison function in (ED above. Assume 
that full regularity and the existence of an uniform upper bound is proved 
for a family of comparison functions w^ a ’ Xa \ 1 < i < n for all r = r(t) for 
(t s ,x s ) £ S ta := {(t,x) £ S\t = t s , } (a section at t s of the CKN-Haussdorff set 
of possible singularities), such that for some m > 2 

V(t 8 ,x a ) £ S ta : w ta ' Xa £ C([r(t 1 ),r(t s )\,H m nC m ), 

(48) 

SU P (t a ,x s )es u su Pr> tl \wI s,Xs \t, .)\ H m <C <00. 

In the following we drop reference to a specific singularity at ( t s ,x s ) for sim¬ 
plicity of notation, i.e., we write 

wi = w t i ‘” Xs , if the reference to (t s ,x s ) £ S is known form the context. (49) 


Here, by an uniform regular upper bound we mean that C is independent of the 
singular point at time t s , i.e., independent of (t s ,x s ) £ St a . For the regularity 
transfer the Leray projection term is crucial. We have 


P,i =P,i 


, dzj 
dxi 


1 


= P,. 


(ts - t ) 


Hence, regularity of Wi, 1 < i < n and (l48ll implies that 


(50) 






(s, x) 


< 


C 


(ts - S)P ' 


(51) 


Recall that we consider p = 1 such that we have the Gaussian G u in ED- a 
main idea for the regularity transfer is to use a spatial symmetry of first order 
spatial derivatives of the Gaussian in convolutions with the Leray projection 
term, where we use this symmetry along with (local) spatial Lipschitz continuity 
of the Leray projection term, i.e., the relation 


At" ^2 (Vm,jVj,m) (s,X-y)-L\„ ^2 ( v m,jVj,m) (s,X-y') 
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< i\y-y I' 

(52) 








for a Lipschitz constant which is independent of x (inherited from spatial Lips- 
chitz continuity of the Leray projection term L\ n ^j m =i { w m,jWj,rrS) ( s > x )i 1 < 
i < n proved in detail below). 

Remark 1.6. Note that local spatial Lipschitz continuity of the Leray projec¬ 
tion term is needed for our purposes, because we estimate convolutions of the 
Gaussian and (mainly of first order) spatial derivatives of the Gaussian. 

We shall prove below, that a refinement of this idea leads to the conclusion 
of a full transfer of regularity from the comparison function Wi , 1 < i < n to 
the original velocity function vt, 1 < i < n. We start with the estimate of 
the velocity components themselves and consider the proof of theorem 11,21 first. 
Furthermore we first construct regular extensions of the function Wi,l < i < n 
to the CKN-Hausdorff set of possible singularities. We then transfer the results 
to the original velocity component functions. 

We get 

Theorem 1.7. Let T > 0 be a given horizon, and let Vi, 1 < i < n be a weak 
Leray-Hopf solution of the Navier Stokes equation, where Vi £ L 2 ([0,T],H l )for 
all 1 < i < n. Then after any finite time the velocity component functions 
Vi,l < i < n and their spatial derivatives up to second order can be continuously 
extended to the Haussdorjf set of possible singularities predicted by CKN-theory. 
The result can be sharpened such that it also holds if a weak Leray-Hopf solution 
Vi, 1 < i < n satisfies Vi £ L 2 ([0,T], H 1 ~ e ) for all t > 0 and e > 0 small. 

The result in (ED does not tell us how the uniform upper bounds depend 
on t s as t s 4- 0. It seems that the analysis of this paper can be extended in order 
to show that H 1 - data (maybe even Lf 1_e -data for small e > 0) are sufficient for 
global existence, uniqueness and smoothness. In any case, local time existence, 
which is available for strong data implies the existence and regularity for all 
time. It is therefore quite possible that a stronger conclusion (with with H 1 or 
H l ~ e instead of H 2 DC 2 data or H^ +1 - data) of CKN-theory than the following 
can be obtained. Nevertheless, we note 

Corollary 1.8. For Vi £ H 1 for all 1 < i < n with respect to time and space 
global existence and uniqueness holds for all time. Moreover if a local time 
contraction results for the initial data, then global existence and uniqueness hold 
for all time t £ [0, oo). This is true for Vi( 0,.) £ H m fl C m ,m >2 on the 
whole space and for i>j(0,.) £ H 2 5 for the torus in dimension n = 3 (sufficient 
criteria, cf. also appendix). 

The local contraction results were proved elsewhere. These local contrac¬ 
tion results are constructive, and they may be of independent value for nu¬ 
merical analysis. The CKN-theory is an example of the additional power of 
non-constructive analysis. Although non-constructive methods cannot be jus¬ 
tified from the constructive point of view, even from this constructive point of 
view non-constructive methods make predictions which are then proved later by 
constructive methods with much more effort (cf. Hilbert’s early paper on polyno¬ 
mial invariants and Konig’s later constructive argument). It is remarkable that 
the statements of Theorem Q and of Theorem H77I are true under the weaker 
condition that a weak Leray-Hopf solution which satisfies Vi £ L 2 ([0, T], 7J 1_e ) 
for any small e > 0. This shows the full power of the CKN-theory combined 
with the singularity analysis, which we consider next in the following proofs. 
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2 Proof of theorem 11.2 


For simplicity we consider p = 1 such that P 2 is a constant and we have G„/ = G(j 
with adjusted viscosity v' = vp. 2 - As a consequence of (1371) we have 

|ui(t s ,x s )| < sup rfoo r \wi (r, 2 j)| 

1 1 I ’ S I I 


< sup Ttoo>ZJ Wi(tin,y)G v '(T,zl - y)dy\ 

+| (E"=i (s,y)G 1 ,/(T-s,2j-j/)d2/ds| 

+ | pi(s)L\ n 1 (wm,jW j>m )) (s, y)G v -(r - s,z T s - y)dyds 


+ 1 Jo fds z (t °- x ‘> w f Z ( s ^y) G ^( T ^ z T s -,s,y)dyds\). 

*1 / 

(53) 

Shifting spatial coordinates we may assume that x 8 = 0 if this is convenient. 
Next we estimate the terms on the right side of (Hull) . The Leray projection 
term is crucial, and we start with this term first. In a first step we consider the 
Leray projection operator more closely. As remarked above it is determined by 
the Poisson equation in (1351) which is the same as the corresponding equation 
for the original velocity component vy by a property of the transformation. We 
prove 

Lemma 2.1. Assume that n < 3 and that for all t £ [ti,t s ] a Leray-Hopf 
solution satisfies vfit ,.) € H 1 . Then 


P,i G L 2 . (54) 

The result still holds under the weaker assumption that for all t £ [ti,t s ] we have 
vfit ,.) £ lL 1_e ° for small eo > 0. Here we note that the upper bound constants 
used are global, i.e., within a fixed arbitrary time horizon T > 0 they do not 
depend on some argument ( t s ,x s ) of a possible singularity. 


Proof. For a Hopf-Leray solution Vi £ L 2 ([t%,t s ], H 1 ') we have 


for all (fix) £ Kf 3 ' Xs) 


\vi,k(t,x)\ < 


(t s 


c 

t)' 5 °|a; — az s | A ° 


(55) 


for some finite constant c £ R+ and for some parameters S, A, which satisfy (case 
n = 3) the relation 

0 < <5 0 < 0 < A 0 < (56) 

The function Vi, 1 < i < n is defined on the whole torus, such that from the 
representation 

Vi(t,x ) := Vi a exp (2iriax), (57) 

aGZ" 
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with time-dependent modes Vi a we get 


x ) =: J2 a ei. » Pai exp (2niax) 

v-' r> - -i Xh fc=l Xl-yez 71 ^ 'yj( a k—7k}vj'yVk(c X — 7 ) /r, . \ 

= Laez- 27ria i l {a7 4 0} ^ h - E” = i^ 2 Qf -- ex P ( 27nQ!X ): 


( 58 ) 


where p^ a denotes the ct-mode of p Note that the infinite vector of time- 
dependent modes vf = (v ia (i))aez ra of the velocity component Vi(t ,.) is in the 
dual Sobolev space of order s £ R iff 


T. I^ia | 2 ( a) 2s < 00 , (59) 

aeZ" 

where 

<a):=(l + |a| 2 ) 1/a . (60) 

Since Vi(t ,.) £ H 1 ,1 < i < n, we have 

y \P\ 2 (v ifs {t)) 2 < oo. (61) 

/3ez« 


Hence, 

47 r 2 7 j(a fe - 7 fc )u i 7 'y fc ( a _ 7 ) < sup c*(t(r)) =: c < oo, (62) 

7 g Z n i 6 [tl,t s ] 

where we may choose (use ab < i (a 2 + 6 2 )) such that 

c*(t(r)) := ^ u 2 7 (t). (63) 

7GZ n 


It follows that 

|p,ia(i)l — £ a gZ" l 27 T* a i|l{a^0} J]"_ 1 47r 2 Q a > (64) 

where p,,;„ denotes the a-mode of p t i- For given t 6 [<i, t s ] the square |p,i a (i )| 2 
has an integrable upper bound, where we note that the factor at in the numera¬ 
tor occurs only on one dimension. Hence for n < 3 we have py(t,.) £ L 2 . Hence, 
since ( t s ,x s ) is the only singular point of p 7 on the cone K l \ tB ’ Xs ' > we have the 
spatial upper bound 


for all (t, x) £ xf s ’ Xs ^ \p yi (t,x)\ < 


— I.T„ — 


(* s - i) 


:|3/2- 


(65) 


where <5i £ (0,1) can be close to 1 because do £ (0,0.5). Finally, if for all 
t £ [ii,i s ] we have Vi(t, .) £ for small eo > 0, then 

\P,ia(t)\ < X/aeZ n l 27r * a ®l l{a/0} 4-7r 2 a 2(1_e,) ’ (66) 


for small e' > 0 such that we still have p,i(t ,.) £ L 2 for n < 3 and the result 


= £ 


a«GZ 


still holds. Here we note again that the factor J2a-e 
J 66 } is finite and the remaining n — 1 dimensional sum is also finite for n < 3 
(comparison with integral upper bounds). 

□ 
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Next, we have the pointwise relation 


H z ■ 

w w 

v ' i ~ v ' i dx i ~ p ' i ( t s -ty ' 


Note that for p = 1 we have 
t 


(tg - t ) 


(tg - t)r = t tgT = t(l + t) t = Y 




+ r 


from which 


tg~t = 


tg(l+T) tg 


1 + r 1 + r 1 + t 


follows. Hence, the function 


xK”3 (r,z) -> \p W i(T,z)\(tg - t) 1 = |p’"(r,z)|^-^ G L 2 


and, since p w is smooth, we have for some 82 G (0,0.5) 


|^M)| < 


(l + r) 1 + 4 a(l + |*| 1 - 5 - e )’ 


(67) 

( 68 ) 

(69) 

(70) 

(71) 


where p™ (r, z) = L\ n (tu m ,j w,-, m )) (r, z). Next recall that for (5 G (0,1) 

we have the standard estimate 


\G v ’{(T,y)\ < 


S |y |n —28 5 


(72) 


or, alternatively, since the integral with respect time is for a > ti n > 0 we can 
even use the obvious estimate 


I G v > (p,y) I 


< 


1 


-3 • 


(73) 


Anyway, recall that Leray projection term is the convolution of the latter two 
terms times \i\ = (t s — t)/t s . Choosing x s = 0 = z T s for r G (0, 00 ) w.l.o.g. we 
get 

SUPrtoo.zJ | It in Jjj / |<iMl( i ( cr ))- £ 'T’* (Ej,m=l K,i®J>)) 0> 2/) X 
x G U ' (r — (7, zj — y)dyda | 

A SUp rtoo z r | f tin J’l y |<! ( 1+(T )2+4i (1^11.5-.) Ir-aPPj-yl 3 - 2 ^^! 


(74) 


< c + sup Tfoo | fi tn f ly]£1 (1+r)1+ ^ 1+ ^ Nj|1 . 5 _ e _ 2 , ) dydgI <c' <oo 

for a finite constant c', and where we may choose <5 G (0,1). Here, we may use 
the elliptic integral estimate 


dy 


Ib » k-y| Q M 6 


< max {c|y| ra-0-b , c} 


(75) 
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for some finite constants c, c > 0 and where is a ball of finite radius around 
x. Such an upper nound of an elliptic integral may be obatined by spliiting the 
integral or by partial integration. 

Similar (simpler) estimates hold for the Burgers term and for the initial data 
convolution term, i.e., we have 

sup r!t00)Zl | ff f ]y |<! (E"=i Mi (wj + Zj) w i,j) (s, y)G v > (r - s,z T a - y)dyds\ < c", 

(76) 

and 

supfl/ Wi{t in ,y)Gv'{r,z T s - y)dy\) < d" (77) 

T,zJ ' Jyefl J 

for some finite constants c", c> 0 by similar considerations as in the case of 
the Leray projection term. Finally, we consider 

sup j / f wf (s,y)G„'(T,zl;s,y)dyds\). (78) 

rtoo.zl Jti n Jd s z) ts ’ Xb ^ ' 

r l 

Here recall that Q s zf ' x “' 1 denotes the spatial boundary of the cylinder zf’ Xs ^ 
with spatial basis fi. Recall that we may assume that 0 = x s = zj for all 
r £ (0, oo). For all points (r, z) on the boundary we have 

wf (r, z) = -2wi\ gSz (. ts ,x a )(T , z) + 2 J n Wi(t ln ,y)G„'(T,z;ti n ,y)dy 
*1 


2(N l * Gj,)(t, z) + ^fc =1 fo In ( %Wi\ dSz (t s ,* s ) (a, £)+ (79) 

\ ti 

2 f n Wi(ti n , y)G v '(a, £; t in ,y)dy - 2 (IV* * G w /)(cr, £)) (t, 2, a, £) rf £dcr, 

where we note that 


G k ,=Gf, and Vi\ SK (t,x) = Wi(T,z)\ gSz ( ts , X3 ). (80) 

Here recall that v l | g denotes the restriction of the velocity component to the 
spatial boundary S'# of the cone Kff s \ According to Hopf’s result we have 
Vi £ L°° ([ii,i s ],L 2 (T”)) , such that for all (t,x) £ we have 


Ui| SK (t,x) 


< 



(81) 


for some 0 < A < 1.5 and some finite constant c > 0. Shifting spatial coordinates 
if necessary, we may assume x s = 0. We may write the boundary of the cone 

,%s) 

S K ~ {(t,x)\t £ (h,t s ) & \x - x s | = t a -1}, (82) 

with t\ £ Uj. This boundary may be written in polar coordinates in order 
to obtain a simple description of the spatial boundary of the corresponding 
cylinder. For ro := t s — t,\ and x s = 0 and p = 1 we have \x — x s \ = \x\ = 
(t s — t) p \z\ = (t s — t)\z\ such that the boundary of the corresponding cylinder is 
described by 

gS z (ts,x s ) |( T; j;)||^| =r 0 & r > tin}. (83) 
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We have x = (t s — t) p z with p = 1, and 


\Gv' (&') y) \qs z^s^s) — 


yf±k 


exp 


(-ML) 

y Ai/'a J 


gSz ^t s ,x s ) 


'/Am 


■ exp - 


( hoi 2 A 

Y Ah''(T J 


Hence, using ifHTTl) . and (IMll . with n = 3 and x s = 0 we have 


(84) 


SU Pr>t in 

— fd s Zi ta ’ 
Z 1 


i(r - cr,0 - y)\ gSz ( ta ,x s ) * G v (cr, y) 

( hoi 2 ^ 
y Ai/'cr J 


ko| A (ts-t(o-)) A 


y/4^ 


exp 


dy 


< d < oo 


(85) 


for some finite c, c*, c' and where we may choose 5 G (0.5,1). The higher order 
boundary terms can be estimated similarly. In the case Vi G H 1 ~ t , 1 < i < n 
for small e > 0 we have observed that Lemma [5] still holds. Furthermore, all 
estimates go through straightforwardly. Finally, note that Vi(ti, .) G C°° (T n ) 
such that the upper bounds can be constructed independently of x s , he., the 
spatial component of a possible singularity (t s , x s ) and the related cone K^ s,x "\ 


3 Proof of theorem 11.71 

First we describe the proof plan. The result that there exists a left continuous 
extension of a Hopf-Leray solution can be sharpened by consideration of spatial 
derivatives using spatial symmetry of first order spatial derivatives of the Gaus¬ 
sian together with local Lipschitz continuity of the Leray projection term. This 
gives the estimates for the crucial Leray projection term convoluted with the 
first order derivative of the Gaussian. The convoluted initial value term and the 
convoluted Burgers term have similar estimates a fortiori. In a second step we 
consider upper bound estimates for spatial derivatives of the boundary terms. 
This argument (given below in detail) implies that for any possible time t s of a 
singularity (t s ,x s ) G S (recall that S is the set of possible singularities predicted 
by the CKN-theory) we have (for some finite constant C ) 

sup \vi(t, .)| H 2 nC 2 < C < °°, (86) 

where [G, t s ) C Uj, and Uj is an open interval where a given Leray Hopf solution 
has full regularity. Furthermore, within a fixed time interval [0,T] for some 
arbitrary time horizon T > 0 the upper bound constants used do not depend on 
the specific location of singularities at (t s ,x s ), i.e., for any sequence (t(?) such 
that (tg,Xg) G S , T* J. t a as k t oo, the finite upper bound constants C t k with 
\ v i{tg , •) \ H 2 nC 2 < G t k have a common upper bound such that sup^^ C t k < 
C < oo. 

First we prove the existence of a regular spatial left-continuous extension. 
The latter task is obtained by proving a regular left-continuous extension for the 
comparison function Wi, 1 < i < n, together with a proof of regularity transfer 
from the function Wi, 1 < i < n to Uj, 1 < * < n. 

We first consider this regularity transfer from Wi , 1 < i < n to i>j, 1 < 
i < n assuming full regularity of the former function and Lipschitz continuity 
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of the gradient of transformed pressure function pf, 1 < i < n, i.e., Lipschitz 
continuity of the transformed Leray projection term. This regularity transfer is 
not trivial (cf. the property of the transformation in (llijl) above). Note that for 
all 1 < j < n we have 

\vi,j(t,.)\ < \wi tj (T,.)\ _ , (87) 

where we again choose p = 1 in the following for simplicity, i.e., in order to 
have solution representations in term of convolutions with standard Gaussians 
at hand. For the regularity transfer the Leray projection term is crucial, and 
we consider this term next. Since r = -d— . or t = f sT A-, we have 

is— t ’ S l+T ’ 


/7 v ■ 

w w X 

’’ ’* dxi ■’ (t s - t) 


T 1 

1 vj / 1 i \ 

= P,i-[ =Pi 


( 88 ) 


We shall observe below that |p“(r, .)| < (even a stronger decay holds). 
Anyway, even by the first two relations in (1881) we have that regularity of Wi , 1 < 
i < n and (1151) implies that 



E i 


u m,j u j,m) 


(s,x) 


< 


C 




(89) 


which is not integrable for p > 1. Nevertheless, we consider p = 1 such that we 
have the Gaussian G% = Gj, in (EU). Note that for p = 1 we have G1 = G v ’ 
with // = j-. where we use the assumption that t s > 0. The function Vi, 1 < 
i <n is defined on the whole torus, say on [—0.5,0.5] n with periodic boundary 
conditions, but we can treat it formally as an initial-boundary value problem 
with artificial boundaries 


d s ([ti,f s ] x T") := {(t,x)\t € [ti,t s ] & Xi = —0.5 or Xi = 0.5, 1 < i < n} . 

(90) 

Let us explain why the Leray projection term is crucial for the regularity transfer 
from the function Wi, 1 < i < n. For the first order spatial derivatives we have 
the representation 

Vi,j(t a ,x) = v i (t 1 ,y)G v>j (t a -h,x- y)dy 

^ S[ tl t s ] xi" (Ej=i(^'Wt) (s,y)G„,j(t a -s,x- y)dyds 

(91) 

"b /[t 1 ,t s ]xT n ■^' T ” i ( v m,ivi,m )^ ( s j y)G v p(t s — s, x — y)dyds 

+ fos atuts] xt -) v f ( s > v)°«d ( T > *; s > y) d v ds - 

Here, the boundary term vf (s, y) is given by an analogous formula as wf 
before. Furthermore, on the right side of m concerning the leading terms only 
the Burgers term and the Leray projection term involve first order derivatives 
of the velocity components Vi , 1 < * < n such that we have to deal with the 
regularity loss expressed in (1871) when we pass form Wij to Vij. This is different 
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for the other terms, where we may use Vi(t, x) = Wi(r, z), and this is true for the 
higher order terms of the expansion of v? as well, of course. The regularity for 
the latter higher boundary terms follows form the regularity of the boundary 
terms for the transformed functions w° straightforwardly. Note here that the 
nonlinear terms (the terms abbreviated by N l in the representation for wf in 
(EHD and (1791) 1 are convoluted twice, and can estimated straightforwardly for 
/z = 1 if a suitable estimate for the Leray projection term is at hand. Hence, 
the Leray projection term is indeed crucial. First note that the assumption 
of a full regularity of the comparison functions Wi on a time interval [fi,t s ] 
implies that we have spatial Lipschitz continuity of the Leray projection term, 
i.e., we obtain this spatial Lipschitz continuity from spatial Lipschitz continuity 
of p™, 1 < i < n, cf. below. This means that for all x we have the relation 


T 1 

-On 


Y ("r. 


.iVj 


(s, x y) Lj 


(v n 


(s,x-y')\ < l\y—y\', 


\j,m= 1 


ij,m=l 


(92) 

for a Lipschitz constant which is independent of x, and for all s£ [t \ 1 t s \ and 1 < 
i < n (inherited from spatial Lipschitz continuity of the Leray projection term 

L{wm,jWj,m)^ (s,x), 1 < i < n proved in detail below). We use this 
Lipschitz continuity together with the symmetry of the first order derivative of 
the Gaussian. In this context for y = (j/i, • • • ,y n ) let 

y~ J = {Vi 3 , ••• yj 3 = -Vj, Vk 3 = Vk for k^j. (93) 


We get (for n = 3) 


ft! /[-0.5,0.5] 3 



(t s - s,. 


y)G v ,j(s,y) 


dyds 


< ft: f- 


ti J[-0.5,0.5] 3 ,y 3 >0 


(■^T" ( V m,j v j,rrSj {t s s j • y) 


-L 


T n 


(E",m =1 ( v m,3 v 3,™)) tys ~ S,X-y J )) 


V.i 


2 (t s —s)\/ A-nvs 




- ex Pl-47|fc) 


dyds 


± fu fi¬ 


ts. J[-0.5,0.5] 3 ,y J >0 


C l\y \ 2 ( |y| 2 

(t.-eWA 6XP (rirfcj. 


dyds 


< ( t fL t )i , for 8 G (0,1), 


(94) 

and for some finite constants C, C *, and where we use the standard pointwise 
estimate 


G l,j( t ~ s ,y) I ~ \u^ G l( t ~ s ,2/)| = 


_ M CXD ( _iKl!—^ 

(t s — s)y/ 47ri/(£—s) \ 4^(t s s) J 


< 


|2 \ <5—n/2—1 ( \ y \2 ^n/2+1-5 


^ (4^(t.- fl) )* lyl (NT n/ " 1 ( 4 J(f- s) ) 


ex P (”47^)) 


< 


— (Aitv(t B -s)) 6 \y\ ri + 1 -' 2b ' 


(95) 
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Note that the latter upper bound is only integrable for S £ (0.5,1), and it is the 
spatial Lipschitz continuity of the Leray projection term which allows to get the 
upper bound in (1M1) with 8 £ (0,1). It follows that for all t £ [ti,t s ] 


/7 

Jti J [—0.5,0.5] 3 



i) (• ^ 8 , 



G v j(s, y)dyds £ L p 


(96) 


for all p > 0 (note that have passed the barrier p = i). Iteration of this 
argument (and application to spatial derivatives) leads to full regularity transfer. 
Hence it remains to prove that we have a) full regularity of the comparison 
function Wi, 1 < i < n, and b) that spatial Lipschitz continuity transfers from 
Wi, 1 < i < n to Vi, 1 < i < n. For the first order derivatives Wi j we have a the 
classical representation 

z) = MUn,y)G p j (T,z- y)dy 

- f z (t s ,*s) (E"= i(^ w o + z i) w i,j) ( s > v)<%,j( T ~ s ’ 2 “ y) d V ds 

(97) 

+ f z (tsfs) yi{s)L\ n (E " m =1 ( w (^>2/)^ i/jjir s^z y)dyds 

+ Jt in f d s z (*s.*s) Wi Z (s,y)GZj(T,z;s,y)dyds, 

tl ,T 

(t X ) 

where Z t v s denotes the cylinder cut off at r > tj. Hence, for all 1 < j < n 

SU Pr to o,ZJ \wij(T,zl)\ 


^ su p (l /{„i ( tl , V )ez^} Wi (***> y) G vA T - * - y)M 

+ | f z (u,x.) (E"=1 i ^ w 3 + Zj)wi,j) (s,y)G P d (T - s,z - y)dyds\ (98) 

+ | fH{s)L l Tn (£" m= i (w mi iwi >m )^ {s, y)G^ j (r — s,z — y)dyds\ 


+ 1 St in fas z (t*,*s) wf(s,y)G^(T, zj; s,y)dyds\). 

We estimate the Leray projection term and remark that the the boundary term 
can be estimated using the ideas of the previous section. Similar estimates hold 
also for the Burgers term and the initial value term a fortiori. Consider the 
cone K ds,Xa associated to a possible singularity (t s ,x s ) £ S and consider a 
Leray-Hopf solution Vi, 1 < i < n on this cone. If Vi £ H l then 


T l 


<j,m =1 


u m,j u j,m) 


(t,x) 


< 


C 


(t s -1) 6 


(99) 


for some finite constant C > 0 and some 8 £ [0,0.5) and small e > 0. If 
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Vi G H 1 e for some small e', then 


y um ’3 u 3> m 

<j,m= 1 


{'^rn,.j ^j.rn) I (f, 


< 


c 


(t s - t) s \x s - x\z 




( 100 ) 


for some finite constant C > 0 and some <5 £ [0,0.5 + e] and small e > 0. We 
estimate the Leray projection term of the transformed equation in the latter 
case. Shifting spatial coordinates we may assume that x s = 0 and have 


(v„ 


(t,x) 


l j,m=l 


< 


c 


{ts - t) S \x\z 


o+e 


( 101 ) 


The left side denotes the first spatial derivative of the pressure with respect to 
the argument Xi . Recall that (with /r = 1) 

P,i = P w i ^ t \ t y or P w i = (t s - t)p ti . (102) 

Recall that for the first order spatial derivatives of the Gaussian we have the 
standard estimate 

\ G ^y)\ < |y|n+l —2^ ' ( 103 ) 

Note again that, alternatively, since <i > tj n . we have a simple estimate 

\Gl tJ (a,y)\<-^-, (104) 

\y\ 

where we may use ^-Gl(a,y) = ^-G\ v (a,y)G\ v (a,y) < C\G\ v (a,y)\ for some 
finite constant C > 0. Hence even stronger estimates than the following hold. 
Recall that Leray projection term is the convolution of the latter two terms 
times y ,i = (t s — t)/t s . A similar reasoning as in the last section leads (in case 
n = 3) to the upper bound 

SU Prtoo I ft, ■/|j,|<l/H( f (°')) I 'T» (Ei,m= 1 (cr, J/) X 

xGl j(r - <7, z - y)dydcr\ 

(105) 

— SU Prtoo | f tl f\y\<l (1+<t) 2 + 4 i(|i/|1- 5 +') | t-o\ b \ t-y^- 26 C ^ 0 ’ I 
— SU Prtoo | ft, (l + T )l + '>l+i C (| z |2.5 + e -25) d<j\ < |.| 0 '.5 + e' • 

for some appropriate constant cq and small c’ > 0. This means that we have for 
dimension n > 2 

SUPrtoo I ff /|y|<l/ i l(*( <T )) i T» (<b y) X 

(106) 

xG i,j( r - y)dyda\ £ L 3_e * 
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for small e* > 0. As the Burgers term and the initial value term have stronger 
regularity, and assuming that the boundary term has the same regularity at 
least (which is indeed true and checked below) we conclude that 

w i°j ( r > •) : = w i,j ( T > ■) G L 3 for r > t in . (107) 

We may use this information as an input in the estimate in (11081) above, i.e., we 
may set up an iterative scheme for k > 1 of the form 

sup T > ti „ K,/(t,.)| 


< sup. 


>ti- 


y\(ti,y)ez' t 


(*s 


s) } 


w[ k 1} {tin, y)G^j(T, z — y)dy 


+ | (E"=i(Mi w f _1) + x j) w ij 1] ) (s,y)G^j(T - s ,, -y)dyds\ 

+| yi{s)L\ n (E" m =i ( s >y) G E( r-s ’■ - y) d v ds I 

+1 ft in w i Xk ~ 1) { s ^y) G u.j{ T ^^ s ^y) d y ds I), 

tl / 


where for the boundary term we have 


a z ,(o) a z 

W, = Wi , 


(108) 

(109) 


and such that w'- is defined recursively and analogously. The functions 
which define the boundary terms inherit regularity vj d ' known from the pre¬ 
vious step and it is easy to check that the upper bounds which hold for the 
Leray projection term are a fortiori upper bounds for the boundary terms (you 
may even use the fact that we are on the boundary of a cylinder where the basis 
is a ball of positive radius around x s =0). We shall observe that gain spatial 
regularity of one order at least at each iteration step. Some embedding results 
may be used here (cf. gj). We have 


g £ H s ' p iff A s g e L p , (110) 

where 

T(A s 5 ) = (1+£ 2 ) s/2 ? (<?)(£), (ill) 

and 

H r ’ q C H s ’ p iff - - - = Eli. (112) 

q p n 

As usual and as before we drop the second superscript in case of L 2 -theory. For 
n = 3 we start with S L 3 and get Wjj(r,.) € H 0 5 . The functions 

Wi, 1 < i < n are defined on a cylinder (not on a torus), but we can adopt 
dual Sobolev spaces obviously (the formal definition may be supplemented by 
the reader). Recall that 


Wi(r, .) G H s iff ^2 \wia{r)\ 2 (a) 2s < oo, (113) 
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If w^ a , a £ Z n denote the a-modes of uj-°- (t, .) then 


K ( °i(r,.)| < -4+7, and |4,a0>-)| < t 4+7 
(a) (a) 


( 114 ) 


for some small e > 0, where we recall 


(a) := (1 + M 2 ) 


2\l/2 


(115) 


Furthermore, the Poisson elimination equation for the pressure is the same for 
Wi as in original coordinates, such that we have 


P W /°\t, z ) =p w i (r, z ) = E Qe z- Pa,i exp (2niaz) 


X)aez n 27r*ail{ 0!f to}' 


M ).;.■= i Id 


lgz „ 4TT 2 'Y J (a k -^ lc )wj~ l w kia 

E?=i 4 ’ r2 «i 


-7) 


exp (27riaz). 


(116) 


It follows that 

4’ (0) 4-)=^(r,.)eH 1 - (117) 

for small e > 0. Iterating this argument it follows that for all k > 1 

.) := Wij{r ,.) e i7 fc_e for r > f in . (118) 

for small e > 0. Hence we have proved a). For b), i.e., the decay of p™ for large 
r of at least order 1 we start with 


n w — T l 

P.i ~ 


J2 Wl ’ 




(119) 


and go back to the representation of the first derivatives of Wi in ©• We use 
again an ’iterative scheme’ based on an equivalent representation. Here iterative 
scheme means that we have indeed a fixed point scheme but iterate regularity 
considerations with respect to this fixed point according to the scheme. We have 

= w i:j (T,z ) = // y | (tjn) „ )ez cw‘s)} wl k ~ 1] {t in , z - y)G^ J (T-t in ,y)dy 
- f z (ts ; * s ) (E"=i(^i w f _1) + x j) w ij~ 1) ) ( t-s,z - y)G^ j (s, y)dyds 
+ / z (ts,*s) in{s)L\ n (E"m=i ( T-s,z - y)G^ j (s,y)dyds 


+ IL f d s z £‘’*•) w i ’ (k 1] ( T ~ s > z ~ y)G^As,y)dyds, 


(120) 

where ; ^ = w m ,i- Consider again p = 1. We observe a gain of regularity 
at each iteration step. As observed above for p = 1 and n = 3 we may us for 
s > ti n the estimate I G+, (s, j/) I < + for some c > 0. All moduli of functions 

I Ij I y/S \y\ 

w m\. = w m,i, = w m , and wf = wf have a finite constant upper bound 
C > 0 for all 1 < l, m < n, such that we get an upper bound for | w ^\| < 
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|Wjn | < ^=, and \wf ’^| < for s > ti n after one iteration step (note that 

the terms with coefficients /xo ~ ~ pi have even a stronger decay by one 

order at this first iteration step. After the second iteration we have the desired 
decay with respect to r. This proves b). Finally we remark that the constants 
inherited from CKN-theory are all global in a given domain with arbitrary finite 
time interval [0, T] with arbitrary finite time horizon T > 0 and do not depend 
on (t a ,x s ) G S. Hence, we have indeed a regular extension of the Leray-Hopf 
solution. 


4 Conclusion 

The Cafarelli-Kohn-Nirenberg theory is one of the powerful tools of the last cen¬ 
tury, which allows to derive conclusions about existence, regularity, and singular 
behaviour for weak function spaces for a considerable class of fluid models. It 
has shown that the treatment of the Navier Stokes equation is hard from the 
perspective of weak function spaces. Conclusions can be obtained which go be¬ 
yond statements of global regularity which assume more regular data. We have 
stated only a few in this paper, but the method outlined here my be used in 
order to investigate the asymptotic singular behavior near weak data as t s J. 0 
etc.. In the appendix an alternative short argument is given which reduces the 
global regularity and existence problem to a local time contraction results for 
H m D C m data for m > 2. We think that the incompressible Navier Stokes 
equation problem belongs to a huge class of evolution problems which have a 
global regular or global smooth solution branch. Many equations of this class 
may have singular solutions, but in case of the incompressible Navier Stokes 
equation a global regular solution branch v G C° ([0, T], H m (~l C m ) for arbi¬ 
trary T > 0 is well-known to be unique, i.e., if fij, 1 < i < n is another solution 
of the incompressible Navier Stokes equation, then we have 

\v{t) - v(t)\ 2 L2 < |«(0) - v{ 0)\ 2 l2 exp (c (K s )Il 4 + K S )|L) d s) ( 121 ) 

where C > 0 with p = 8 in dimension n = 3. From the perspective of CKN 
theory after small time to > 0 we can close the gap between 

Vi G L°° ([t 0 , T],L 2 (T”))n L 2 ([0,T],Ft 1 " e (T n ))nL 8/3 ([0,T],L 4 (TP)) small e 

and 


vt G L°° ([t 0 , T],L 2 (T")) n L 2 ([0, T], H 1 (TP)) n L 8 ([0, T\ , L 4 (T n )) , 

and the so-called global regularity and existence problem follows then from local 
time contraction, if strong data, say Vi( 0,.) G H m (~l C m , m > 2 are assumed. 

5 Appendix 1: Comparison to another global 
regularity argument 

We consider the reduction of the global regular existence problem to a local 
time contraction result. This argument is considered for the problem on the 
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whole domain R". It is a variation of arguments given elsewhere, and works 
even without viscosity damping estimates. The argument can be reformulated 
for the n-torus. Local time contraction on a short time interval [to, to + A] for 
data v%(to, .), 1 < * < n, to > 0 with Vi(to, .) £ H m n C m for m > 2 shows 
that there is a local- time representation of the velocity component functions 
v\, 1 < i < n of the form 


< = «<o, ■) «.p Gp - E“.i (»J Jf ) * C* 

+ (e”„„/ r , (&*»<-»>)E”„-1 (S“£) (•.»)<*») *G, 

Using the incompressibility condition 


£^= 0 , 


j=i J 


( 122 ) 


(123) 


we may rewrite the Burgers term, where we have 

3(W) = 

i=i 


E- 


9a;, 


Tl r\ Tl Tl r\ 

E OVi GW 7 a Vi 


9a;,- 


9a; 


a=i J i=i J i=i J a'=i 
Hence the local representation in (11221) may be rewritten in the form 

< = <(to, ■) - £"=i * G Vtj 

+ (/,- (*»(■ - V)) E" m =l (•■»)<*») * 


(124) 


(125) 


where all nonlinear terms are convolutions with a first order spatial derivative of 
the Gaussian G„. We may use the Lipschitz continuity of the Leray projection 
term function for strong data, which we get close to data at to > 0 by local time 
contraction in regular space with respect to the norm sup te [ toto+A ] \v"(t, -)\ HmnC , 
for m > 2. For the first order spatial derivatives of the Gaussian we compute 


Gv,i it, y ) 


-2j/i 1 

47 list V47T!/i n 



< 1 (\jAl\ 5 ~ n/ 2 

— (Anist) 6 \y\ yAms J 


|2 \ n/2+1-5 


exp 



(126) 


Hence we have for S £ (0,1) 

\G v ,i(t,y)\ < 

where the constant 


C 


{Aitvt) s \y\ 


5 | 7 ,|n-(-l —2<5 : 


C = sup (. z ) n/2+l ~ s exp (-z 2 ) > 0 
M>° 


is independent of v > 0. Similarly we get 

\G r v At,y)\ < 


C 


2+\5\„.\n+l-25 ’ 


(47T^r 2 t)^|7/| 


(127) 


(128) 
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for the scaled equation with transformation v\' v (t, y) = v\ ( t , x) and y = rx. For 
1 < \/3\ < m we have 


D&7 = D^(t 0 , .) * sp Gl - r E ; =1 Dl {?YvY) * G r vj 

+ r (/*« ~ y)) Ej,m=l D X (;V)dy) * G r vJ , 


(129) 


where for |/3| > 1 we have 0 < |'y| < m — 1 , 7 ; + 6ij = /3i for all 1 < l < n, 
and where we used the convolution rule. Now consider the comparison function 
uf 0 , 1 < i < n where for all y = rx 

t-t 0 


(! + t)ui ,to (s,y) = Vi{t,y) = Vi(t,x), s = 


(130) 


V 1 - (t - to) 2 

We may consider this transformation for data i>l(io,-) on a time interval t — 
to = 0.5 which corresponds to a time interval of length 0, -^= in terms of the 
transformed time s. Next by local time contraction results we have the local 
representations 

u r i' to {s,x ) = f Rn u r i ' to (Q,y)G l Z’ r (s,x-, 0 ,y)dy 
~ Jo Jr n M(o)<’ to (c, y)G^ r (s, x; a, y)dydo- 

~ fo /r" ry T ’ 2 {a) E"= 1 (uJ’ to <’*°) (cr, y)G£J(s, a:; cr, y)dydcr ^ 3 ^ 

+ fo / R n r/r r ’ 2 (a) E£ r =i / R „ (#«(* - i/)) x 
'^r ’* 0 duZ’ to ' 


E n 

. 7 . 1=1 




(o, J/)G£[(s, a:; cr, z)dydzda , 


and for the multivariate spatial derivatives of order 1 < |/3| < m we have 
D x u J t0 ( s » ®) = Jr n DYE 0 (0, y)G^ r {s, x; 0, y)dy 

- /o’ /r" v( a ) u i?p ( CT > y)GZ' r (s , x; cr, y)dyda 

- fo S /r" r M r, 2 (o) E"=i { u T° u J t0 ) ,0 ( CT ’ 2 /) G £J( s > 

+ fo /r» G^’V) E",r=l /r» CM* - 2 /)) X 

'a«r ’ t0 9<’*° 


(132) 


E n 

. 7 ) 1=1 


dx« dxi 


(o, y)G^’[ (s, x; cr, z)dydzda 1 


where we use the integration of the Burgers term above and we have for s £ 
0 , and for fc £ { 1 , 2 } we have 


y = 


u ( s ) - Vi-(t(«)-«o) 23 > 3V3 _. 

— l+t(s) — 8 ( 1 +T) — • HO, 


r ^T,fc r (q _|_ t(s)) k y = r-y /1 - (i(s) - to) 2 (1 + t(s)) 
< r(l + T). 


(133) 
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(134) 


Here, G^ ,r is fundamental solution of the heat equation 


GZ - /A G^ r = 0. 


Now the local solution representation with first order spatial derivatives of the 
Gaussian has a Levy expansion with leading term upper bound (needed on a 
compact domain) 




C 


(yovr 2 {s — cr)) 5 1 y — z\ n + 1 ~ 2S ’ 


( 135 ) 


where C > 0 is a finite constant (dependent only on dimension). Note that the 
term G^f(s, y; a, z ) is not integrable for S £ (0, L) but convolutions l * G^f with 
a Lipschitz continuous function y —> l x (y) = l(x — y ) are, i.e., for such a function 
with Lipschitz constant Iq there is some finite constant C > 0 (dependent only 
on dimension) such that on a ball B of radius yovr 2 around the origin the 
increment over a time interval [to, s] of such a term has the upper bound 


fto SU P b I G^l 0, y) | < f t s o J B ^ 0Vr * { J°? 0 y\ z \n- 2 > dzdcj 

< l 0 C{y 0 vr 2 ) s (s - to) 1-5 , 


( 136 ) 


where Iq is a Lipschitz constant of l x . In the case of data Vi(to,.) £ H m D 
C m for to > 2 at time to > 0 a local time analysis over a msall intervalk 
[to, to + Ao] shows that the Burgers term functional and the Leray projection 
term functional and their firsr order derivatives have a Lipschitz constant loA-g 
which is independent of x. 

Remark 5.1. For (113611 we may use the observation 


ft 0 su pr- \ l ( x ~ y) G uf O’ V ) I 

= ft 0 SU PR” \Kx-y)j^G^ r {s,y)\ 

< ftl sup R n >y .> 0 \\l x (-y) -Zx(-y _,i )l|^:G')’ r (s,2/)| 

< ftl su Prm, 4 >o \lo^G^ r {s,y)\ 

where y 1 ~ = [y l {~, ■ ■ ■ , yh~) with y-~ = y i for j ± i and y\~ = -y t . 
For a small time interval [io,s] the complementary integral 



( 137 ) 


( 138 ) 


becomes relatively small, i.e.,for every e > 0 and Ao = (s — to) small enough we 
have 


ft S 0 fw n \B \l* G ^J(s,y)\dyds 


< el 0 C(y 0 isr 2 ) 6 (s - to) 1 5 


= el 0 C(y Q vr 2 ) s Al s . 


( 139 ) 


This holds because the exponent of the Gaussian becomes small (note that the 
radius of B does not depend on time - the estimates are designed for the case 
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v > 0). Recall that l is Lipschitz such that that the integrand of (11381) has an 
upper bound 


l*Gv,i(s-to,y, 0,0) <c^rexp( 


-A 


\v \ 2 

4:/j, 0 vr 2 (s—t 0 ) 


)• 


(140) 


where \y \ 2 > povr 2 and A > 0 is a finite positive constant which depends on 
Ho > 0. Hence if 

r < a/Ao °, for e 0 > 0 small (141) 

then for any e > 0 the relations in (11361) and in (11391) holds for a time step 
size Ao > 0 small enough. This implies that solution increment due to the the 
nonlinear terms have an upper bound which is proportional to Aq^ 1 e °l <5+1 s 
where op : = 2(1 — eo )S +1 — S > 1 for e 0 small enough and given 6 £ (0,1). This 
means that the increment Svi := Vi — Vi(to, ■) *sp G v (* sp denoting convolution 
with respect to the spatial variables) has an upper bound 


max sup \Svi(cr, .)\ H m nC m < A a ° (142) 

l<i<n a£[t 0 ,t 0 +A 0 ] 


with some ao > 1 on a small time interval Ao, a growth which is offset by the 
potential damping term over the same time interval if Ao is small enough. 

Next local contraction with respect to a H 2 D C 2 -norm implies that we 
have Lipschitz continuity of the Leray data function and its first order spatial 
derivatives. More precislcly 


= C 2 < 00 there 


Lemma 5.2. Given data u r /°(to, ■) with u r /°(to,-) 

J J H 2 nc 2 

exists a time step size Ao such that u r j to (s ,.) € H 2 D C 2 for s £ [to, to + Ao]. 

Moreover, the function 


+ E",r=i Jr™ ( Kn,i{y - z)) Ujf r °(s, z)u r r ’j(s, z)dz 


(143) 


is in C 1 D H 1 . 


Proof. The first statement follows form local contraction of a standard iteration 
local solution scheme with respect to a H 2 D C 2 norm (supremum over local 
time). Concerning the regularity of the Leray data function we first observe 
that for functions u r f to (s, .) £ C 2 D H 2 such that the first spatial derivatives are 
in C 1 and globally bounded. Hence we have 


-\ 7 > 


{ s ,.)u r r f;(s,.)£C 1 nL 2 cHl 


(144) 


where Hf oc denotes the Sobolev space which is locally H s for exponent s £ 
R (Sobolev L 2 -theory). As data are in Hl oc we know from the regularity of 
uniformly elliptic operators of second order that u\ o (s ,.) £ = Hf oc . Hence, 

{s,.) £ Hf oc . For dimension n = 2 it follows that the right side is indeed 
in H 2 oc (product rule for Sobolev spaces. Hence for data in H 2 D C 2 we have 
indeed a regular Leray data function such that the first order spatial derivatives 
are Lipschitz. 

□ 
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( 145 ) 


Hence given a time horizon T > 0 for the choice 

1 


1 + T 


and Lipschitz continuity (Lipschitz constant Iq) of the (first order spatial deriva¬ 
tives of)Euler-Leray data function, i.e., Burgers term operator plus Leray pro¬ 
jection form operator applied to data u[’ to (s,.), the nonlinear terms in (113211 for 
n > 3 have an upper bound 


Crn(s)(4'K^(s)i / r 2 ) < 


1 


1 + T 


1+25 


, Se (o,i). 


(146) 


Note that the damping term (i.e. the potential term, or the second term on the 
right side of (11321) 1 has no parameter r and the damping over a time interval Ao 
is of order ^Ao times the data norm (for small Ao). More precisely, we have 

Lemma 5.3. Given a time horizon T > 0 and |it[’ to (0, -)| ffmnCm < C m there 
exists a constant c(n , m) > 0 depending only on the dimension n and the regu¬ 
larity order m > 2 and a parameter 

c(n,m)(C m + l) 2 (l + T) 

such that the representation in \1S2\) holds on the time interval 
0 < |/3| < m and we have for all s £ 


73 


for 


0 ’ V3 


r,t 0 


(s, ■) 


H mn C m < K 


r,t 0 


( 0 ,.) 


ff™nC" 


< C„ 


(148) 


Remark 5.4. The constant c(n,m) can be computed explicitly and contains 
upper bounds of local L 1 of K, as factor of one of its summands etc. 

Now let a time horizon T > 0 be given, and assume that for the paramter 
r > 0 of the preceding lemma 


(1 + to)C m > (1 + fo)|'M[’ t °(0, .)| 


_ I r,to 


(*0, -)| 


H™nc n 


(149) 


H m nc m ~ l 

has been proved up to some time to > 0. Then according to the Lemma above 

we have 


for this r > 0 and for all s £ 
K to (s,.y 


0) 73 


H m nc m ^ l M 


r,t 0 


( 0 ,.) 


- Crr 


Hence for t £ [to, to + 0.5] we have 


r,t 0 


(*,0| 


<(l+t)<^(0,.) 


H m nc ri 

r,t 0 . 


< (l+f(s))K’ to (s,.) 


H m nC” 




< (1 + t)C m . 


(150) 


(151) 


Hence, for given time horizon T > 0 there exists r 
0 < t < T 


> 0 such that for all 


I r ,to 

k- 


k-) 


lHmnCm <(l + t)C m . (152) 

We observe that this argument depends on the (mathematically) strong as¬ 
sumption H m D C m , m > 2. The CKN-theory provides this situation for any 
small time, and it tells us about the asymptotic behaviour of singularity upper 
bounds near the T 2 -data. 
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